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ABSTRACT 


V 

A  wide  range  of  passive  optical  sensor  applications 
often  require  the  sensor  to  operate  within  the  atmosphere 
while  the  objects  it  examines  are  outside  the  atmosphere.  The 
refraction  of  light  by  the  earth's  atmosphere  becomes  a 
significant  error  source  when  the  objects  are  at  low  elevation 
angle.  When  the  measurement  accuracy  at  low  elevation  angle 
is  important  for  the  sensor  application,  an  accurate 
refraction  compensation  scheme  is  needed.  In  this  report,  we 
provide  an  algorithm  that  will  compensate  the  error  when  the 
object  is  at  a  finite  range  from  the  sensor.  Sensitivities  of 
compensation  errors  to  object  range  error,  atmospheric  model 
mismatch,  and  possible  statistical  variations  about  a  given 
atmospheric  model  are  also  provided. 
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I.  INTRODUCTION 

The  refraction  of  light  by  the  earth's  atmosphere 
has  been  observed  by  astronomers  for  centuries.  In  order  to 
obtain  the  actual  star  location,  one  needs  to  know  the  path  of 
the  light  ray  from  the  star  through  the  atmosphere.  Mathemat¬ 
ical  formulas  that  trace  a  light  ray  through  nonuniform  medium 
have  been  derived  and  can  be  used  for  this  purpose . ( 1 1  - ( 2] 

With  limited  atmospheric  data  around  the  earth,  it  was  impos¬ 
sible  to  compensate  the  bending  of  light  accurately.  Until 
recently,  numerous  data  have  been  collected  around  the  earth 
in  different  seasons  and  at  different  altitudes.  Atmospheric 
models  have  been  generated  and  used  to  correct  these  errors, 
for  example  six  atmosphere  models  are  provided  in  a  computer 
code  for  atmospheric  transmittence  and  radiance  computation 
called  LOWTRAN  5.(3J 

A  wide  range  of  passive  optical  sensor  applications 
often  require  the  sensor  operate  within  the  atmosphere  while 
the  objects  it  examines  are  outside  the  atmosphere.  When  the 
measurement  accuracy  at  low  elevation  angle  is  important  for 
the  sensor  application,  an  accurate  refraction  compensation 
scheme  is  needed.  There  are  computer  codes,  for  example 
LOWTRAN  5(3],  which  can  be  used  to  calculate  the  refraction 
error  and  compensate  it.  It  is,  however,  not  applicable  to 
the  case  when  an  object  is  at  a  finite  range  from  the  sensor. 
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Figure  1.1  shows  that  for  each  true  target  position,  the 
apparant  elevation  angle  can  be  determined  uniquely.  However, 
given  the  apparent  elevation  angle,  the  true  target  elevation 
is  determined  only  if  the  range  to  the  target  is  known.  This 
is  the  purpose  of  this  report  to  provide  an  algorithm  that 
will  generate  the  actual  refraction  error  for  objects  with 
finite  range  and  also  provide  some  sensitivites  of  these 
errors  to  the  object  range,  the  underlying  atmospheric  model 
and  possible  statistical  variations  about  a  given  atmospheric 
model . 

The  report  is  organized  in  five  sections.  Section  2 
provides  a  brief  review  of  Fermat's  principle  and  derives  an 
important  formula  for  ray  tracing  in  spherical  coordinates. 

In  Section  3,  a  specific  form  of  refractive  index  (a  ratio  of 
velocity  of  light  in  vacuum  and  in  a  medium)  is  assumed, 
namely,  it  is  a  spherically  symmetric  function  about  the  cen¬ 
ter  of  the  earth.  The  formula  that  computes  the  refraction 
error  is  derived  for  an  object  with  a  range  R  from  the  ob¬ 
server.  A  simplified  version  of  Edlen's  expression  for  the 
refractive  index  of  air  [4]  is  used  in  Section  4  to  generate 
the  refractive  index  of  air  at  different  altitudes  from  at¬ 
mospheric  profiles  provided  by  LOWTRAN  5.  The  sensitivities 
of  compensated  elevation  refraction  errors  to  the  object 
range,  atmospheric  model  mismatches  and  the  statistical  varia- 


II. 


RAY  TRACING  IN  SPHERICAL  COORDINATES 


When  light  rays  are  propagated  through  an 
atmospheric  medium  of  continuously  varying  refractive  index, 
they  experience  a  change  in  direction  or  refractive  bending. 

In  this  section,  we  will  derive  the  equations  for  ray  tracing 
in  spherical  coordinates  from  Fermat's  principle. 

2.1  Fermat's  Principle  M 1 

Let  n  be  the  refractive  index  of  a  light  ray  in  a 
medium.  The  optical  lenqth  of  a  ray  which  joins  points  P1  and 
P2  (denoted  by  [P.|P2])  is  given  by 

P2 

[p1p2l  -  Jp*  n  ds  (2.1) 

where  s  is  the  distance  measured  along  the  light  ray. 

Fermat's  principle  states  that  a  light  ray  always 
chooses  a  trajectory  that  minimizes  the  optical  length.  In 
mathematical  terms,  Fermat's  principle  assumes  the  form, 

P2 

[P1P2J  =  / p  n  ds  =  minimum.  (2.2) 

Instead  of  the  optical  length,  we  can  introduce  the  concept  of 
transit  time  by  dividing  (2.2)  by  the  constant  c,  the  velocity 
of  liqht  in  vacuum.  Since  the  velocity  of  light  in  a  medium 
with  refractive  index  n  is  c/n,  we  have 
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(2.3) 


P2 

(P1P2]  =  c  / p‘  dt. 

Therefore,  Fermat's  principle  is  also  known  as  the  principle 
of  least  time. 

2.2  The  Euler  Equations  in  Spherical  Coordinates 

Let  us  consider  that  the  refractive  index  is  a 
smooth  and  continuous  function  of  position  in  spherical 
coordinates  and  is  denoted  by  n(r,6,$).  We  use  the  definition 
of  the  line  element 

ds  =  /dr2  =  /r2de2  +  r2sin29dij>2 


-  /I  +  r20 ' 2  +  r2sin2e<t> '  2 


dr 


with 


fl.  =  li 

6  d? 


i  _ 

dr 


to  express  (2.1)  in  the  form 


/ p  L(r ,  9 , 4> ,  9  ' ,  <(> '  )dr  =  minimum. 


(2.4) 


(2.5) 


(2.6) 


The  function  L  is  given  by 

L(r,9 , 4>  ,0  '  ,<(> ' )  «  n(r,9,4>)/1  +  r20'2  +  r2sin294>'2  (2.7) 

In  variational  calculus,  the  function  L  in  (2.6)  is  called  the 
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Lagrangian.  The  solution  of  (2.6)  is  well  known  and  can  be 
obtained  by  solving  the  Euler  equations: 


d 

dr 

a  l 
a  e ' 

_  Ik  =  o 

34)  Uf 

(2.8) 

d 

dr 

aL 
a  $ ' 

3L 

34i  “  °* 

(2.9) 

Substituting  (2.7) 

into 

(2.8) 

and  ( 2.9)  ,  we 

have 

2 

d  nr  0 

_  rr  c 

,  an 

2  . 

nr  sm0cos0 

,2 

(2.10) 

dr  s ' 

—  o 

30 

+  s' 

2  2 

d  nr  sin  0 *  *  ,  an 

dr  s '  a  <j> 


with 


s 1  =  >  1 

s  dr  -  '• 


(2.11  ) 


(2.12) 

With  (2.12)  and  the  chain  rule,  these  equations  can  be  written 


as 


d  2  de ,  an  .  2  .  n  . .d* . 

d?  (nr  di}  =  aT  +  nr  sinecosetgf) 


(2.13) 


d  2  2„  d <{> .  an 

a?  (nr  Sln  9  a?1  -  ?♦  • 


(2.14) 


These  two  equations  should  be  sufficient  to  determine  the  ray 
trajectory.  The  corresponding  Euler  equation  for  r  can  be 
derived  from  (2.13)  and  (2.14)  using  the  constraint  (2.4). 
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2.3  The  Spherically  Symmetric  Case 

For  the  case  where  the  refractive  index  depends  only 
upon  r,  (2.13)  and  (2.14)  become 

3s  *nr2  Hi-)  =  2sine cos 8  )  2  (2.15) 

^  (nr2sin29  =  0.  (2.16) 

Integrating  (2.16),  we  have 

nr2sin20  =  C^.  (2.17) 


We  may  choose  the  coordinate  system  so  that  ^  *  0  initially. 
Then,  Cj  -  0  and  (2.17)  yields 


for  all 
problem. 


d$ 


ZJ 


*  o 


s.  This  reduces  the  problem  into  a  two  dimensional 
that  is 


a?  (nr2  S?*  “  °*  (2-18) 

This  result  will  be  used  in  the  following  sections  to 
calculate  the  error  of  angular  measurement  due  to  refraction 
in  an  atmosphere  when  spherical  coordinates  are  specified  and 
the  refractive  index  of  the  atmosphere  is  given. 
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III. 


ELEVATION  REFRACTION  ERROR  FOR  SENSOR  VIEWING 
THROUGH  ATMOSPHERE 


In  this  section,  we  will  limit  ourselves  to  describ¬ 
ing  the  gross  characteristics  of  atmospheric  refraction. 

Other  causes  of  refraction  such  as  irregularities  and  varia¬ 
tions  within  the  earth's  atmosphere  which  can  not  be  deter¬ 
mined  by  theory  will  not  be  considered.  Furthermore,  we 
assume  the  refractive  index  of  the  earth's  atmosphere  is 
spherically  symmetric  and  depends  only  upon  r,  the  distance  to 
the  center  of  the  earth.  The  non-sphericity  of  the  earth  can 
be  taken  into  account  locally  while  applying  results  of  this 
section  by  using  the  effective  earth  radius  at  a  given  lati¬ 
tude  and  the  associated  model  atmosphere. 

Let  us  consider  a  light  ray  passing  through  PQ  and 
P1  in  the  earth's  atmosphere  as  shown  in  Figure  3.1.  Due  to 
the  spherically  symmetric  assumption,  the  trajectory  of  the 
light  ray  lies  in  a  two  dimensional  plane  and  can  be  repre¬ 
sented  in  polar  coordinate  as  shown  in  Figure  3.1.  Further¬ 
more,  the  trajectory,  s(r,9),  satisfies  (2.18),  or  equivalent¬ 
ly* 


nr 


2  de 
ds 


constant. 


(3.1  ) 


Let  us  define  the  ray  inclination  angle,  41,  to  be 
the  angle  between  the  tangential  vector  of  the  ray  and  the 
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OPTICAL  PATH 


Fig.  3.1  The  trajectory  of  a  light  ray  lies  in  a  two 
dimensional  plane. 


local  horizon.  Then,  (3.1)  becomes 

n  r  cos  ip  =  Co  (3.2) 

where  CQ  is  determined  from  the  initial  values  of  n,  r  and  \p 

at  Pq,  that  is 

n  r  cos  (p  =  nQ  rQ  cos  ipQ .  (3.3) 

If  the  refractive  index  were  constant  between  Pq  and  P^ ,  (for 
example,  the  atmosphere  removed),  (3.3)  reduces  to 

r^  cos  ip^  =  rp  cos  \pQ.  (3.4) 

The  equality  holds  when  the  ray  is  a  straight  line.  The 
angle,  >Pg(  =  |P^)  is  commonly  called  the  elevation  angle  of  an 
object  at  P1  with  respect  to  an  observer  at  PQ .  The  range,  R, 
is  defined  as  the  distance  between  PQ  and  Pj . 

As  demonstrated  in  Figure  3.2,  an  object  at  P1 
(r=r.j)  with  a  range  R  from  PQ  would  be  observed  to  have  an 
apparent  elevation  angle  <Pq  at  Pq  ,  as  if  it  were  seen  at  P1  * 
with  a  range  R  from  PQ .  If  no  atmosphere  were  present,  the 
true  elevation  angle  would  be  ot  as  indicated  in  Figure  3.2. 

The  elevation  error  due  to  refraction  is  defined  as 

Aa  =  ipQ-a.  (3.5) 

Let  ip^  be  the  angle  between  the  tangential  vector  of 
the  ray  at  P^  and  the  local  horizon  as  shown  in  Figure  3.2. 
Then,  we  have 

a  =  ip^  -  6  (3.6) 
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139730-1 


where 


ip  i  =  cos 


) 


1 


and 


(3.7) 


3 


de  , 


(3.8) 


with  Cq  =  r^  cos  ipg  and  ng,  n^  are  refractive  indexes  at  Tg 
and  r.,  respectively. 

j  „ 

Using  the  fact  that  sin  ip  =  (3.1)  can  be  re¬ 

written  in  the  form 


2  .  de 

nr  simp  gp  =  CQ .  (3.9) 

Substituting  (3.9)  into  (3.8)  and  sinip  by  a  function  of  cosip, 
and  using  (3.2),  we  obtain 


-  /? 


co  dr 


nr 


2fT- 


(Cg/nr) 


(3.10) 


Notice  that  the  refractive  index  n  is  a  function  of  r.  There 
is  no  closed  form  expression  of  this  integral  in  general. 
Whenever  the  function  n  is  specified,  6  can  be  determined  by 
numerical  integration  for  a  given  ipQ. 

For  a  given  set  of  parameters  (a,R,rg),  Aa  is  deter- 


mined  by  solving  <|/g  such  that 


rising  =  R  cosa 


(3.11  ) 


where  e  is  a  function  of  as  given  by  (3.10)  and 


ri  8  V  R  +  rn  +  2Rrrtsina. 


For  the  case  when  n  -  nQ  over  [r^r^,  we  have 


0  «  cos~1(C0/n0r1)  -  cos'1(C0/n0r0) 


(3.12) 


*  “  *1  “  V 


(3.13) 


It  is  easy  to  see  that  Ao»0  and  (3.11)  is  satisfied. 


IV. 


REFRACTIVE  INDEX  OF  AIR 


As  mentioned  in  Section  I,  the  refractive  index,  n, 
of  a  medium  for  a  light  ray  is  defined  as 

n  =  v/c  (4.1) 

where  c  is  the  velocity  of  the  light  propagating  in  vacuum 
while  v  is  the  velocity  in  the  medium.  The  velocity  of  a 
light  ray  propagating  through  the  atmosphere  varies  with 
changes  in  atmospheric  composition,  pressure,  and  tempera¬ 
ture.  It  is  strongly  wavelength  dependent  at  optical  wave¬ 
lengths,  but  it  is  not  affected  appreciably  by  water  vapor. 

In  this  report,  we  adopt  a  simplified  version  of  Edlen's 
expression  for  the  refractive  index  of  airf4]  as  follows: 

n  *  1  +  10~6  x  (77.46  +  0.459/X2)  |  (4.2) 

where  \  =  wavelength  in  micrometers  (ym) 

P  =  atmospheric  pressure  in  millibars 
T  -  atmospheric  temperature  in  deqrees  Kelvin. 

The  atmospheric  temperature  and  pressure  variations 

are  principally  functions  of  altitude,  season  and  latitude. 

4.1  Model  Atmospheres 

Six  model  atmospheres  used  in  LOWTRAN  5^1  (a  com¬ 
puter  code  for  atmospheric  transmittance  and  radiance  computa- 
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Pig.  4.1  The  temperature  profiles  as  a  function  of  altitude 
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tion)  are  considered  in  this  report.  The  model  atmospheres 
correspond  to  the  1962  U.S.  standard  atmosphere  and  five 
supplementary  models;  that  is,  Tropical  (15*N),  Midlatitude 
Summer  (45°N,  July),  Midlatitutde  Winter  (45#N,  January), 
Subarctic  Summer  (60°N,  July  ),  and  Subarctic  Winter  (60°N, 
January).  The  different  models  are  digitized  in  1-km  steps 
from  0  to  25  km,  5-km  steps  from  25  to  50  km,  then  at  70  km 
and  100  km  directly  as  qiven  by  McClatchey  et.al.,15]. 

The  temperature  profiles  for  these  models  as  a 
function  of  altitude  are  shown  in  Figure  4.1.  The  pressure 
profiles  are  given  in  Figure  4.2.  Figure  4.3  shows  the 
profiles  of  pressure  to  temperature  ratio  for  these  models. 
Notice  that  the  1962  U.S.  standard  model  has  a  pressure  to 
temperature  ratio  profile  very  close  to  the  mean  of  these 
profiles.  The  refractive  modulus 

N  =  ( n— 1 )  x  106  (4.3) 

which  is  proportional  to  the  pressure-to-temperature  ratio 
according  to  Equation  (4.2),  is  also  given  in  Figure  4.3  at 
X  =  1 1 gm. 

In  the  next  subsection,  we  will  describe  the  scheme 
of  calculating  the  refraction  error  for  an  object  at  a  given 
range  and  elevation  angle.  Sensitivities  of  refraction  errors 
to  the  object  range,  atmospheric  model  mismatch  will  be 
presented . 
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ALTITUDE  (km) 


Pig.  4.2  The  pressure  profiles  as  a  function  of  altitude. 


ALTITUDE  (km) 


N 


Fig.  4.3  The  profiles  of  pressure  to  temperature  ratio  as  a 
function  of  altitude. 


4.2  Elevation  Refraction  Error  Computation  For  A  Given 
Atmospheric  Model 

In  this  subsection,  we  apply  the  refraction  equations 
and  various  atmospheric  models  to  show  how  well  the  true 
target  elevation  can  be  estimated  from  the  apparent  target 
elevation  and  approximate  knowledge  of  the  target  range  and 
atmospheric  parameters. 

Let  us  assume  that  the  refractive  index  of  a  given 
atmospheric  model  for  a  light  ray  of  wavelength  XQ  can  be  ex¬ 
pressed  by  the  following  piecewise  constant  function 

n(h)  =  n.  for  h.  .<  h  <  h.  (4.4) 

where  h  is  altitude  in  km  and  h^'s  are  altitude  steps  defined 
in  the  LOWTRAN  5  code.  Substituting  (4.4)  into  (3.10),  we 
have 


cos-1  lc0/nM+1 ( Re+hM ) 1 
cos"1 [C0/ni+1 (Re+hi)]J  (4.5) 


20 


p  =  cos"1 (C0/nM+1(Re+hT)]  - 

+  t  {cos"1  [C0/n.+1(Re+hi  +  1)]  - 
l  =  L  ' 

+  cos"1  [CQ/nL(Re+hL)]  -  ij)0 


where 


C0  =  nL(Re+hs)COS,|'0  (4.6) 

R  is  the  earth  radius  and  h  ,  h  are  the  sensor  and  target 

6  SI 

altitudes,  respectively.  The  indices  L  and  M  are  determined 
by  hs  and  h<j.  as  below: 


h 

h 


r  1  <  h 

L- 1  s 


M  -  hT 


h 

hM+1 


(4.7) 


For  a  given  set  of  parameters  (a,R,h  ),  an  iterative  scheme 
has  been  used  to  solve  for  ij^  such  that  (3.11)  is  satisfied, 
that  is 


(Re  +  hT)sinB  =  R  cosa 


(4.8) 


where  3  is  a  function  of  ipg  as  given  in  (4.8)  and  hT  is  ob¬ 
tained  by  the  following  relation: 

hT  =  [(R0+hs)2  +  R2  +  2R(Re+hg)sina] 1/2  -  R0.  (4.9) 

As  illustrated  in  Figure  4.4,  the  iterative  scheme 

used  in  this  report  is  described  as  follows: 

1.  For  a  given  set  of  parameters  (a,R,h  ),  determine 

hT,  L  and  M  according  to  Equations  (4.9)  and  (4.7), 
respectively.  Let  <|ig  =  a  and  j  =  1  . 
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2. 


Compute  3j  by  (4.5)  for  the  given  \j>0. 
A  new  target  range  is  computed  as 


3. 


R.  =  [(R  +h  )2  +  (R  +h_)2  -  2 ( R  +h_ ) ( R  +h_  )cosf?  . ] 1 /2 


5. 


e  s'  e  t  “ '  *'e ’ "s  ' '  "'e  ’  "T  ‘ 

and  the  associated  elevation  angle  as 


-1 


Oj  =  COS  [  (Re+hT)sin3;./Rj  ]  . 


The  difference  of  a  and  aj  is 


A  a  .  =  a  -  a  .  . 
3  3 


4.  If 


A  a  j  _<  e  (tolerance  limit  of  the  algorithm) 


go  to  step  5.  Otherwise  set 

j  =  j  +  1 


*0  =  *0  +  &otj 


(4.10) 


and  go  to  step  2. 

The  elevation  error  due  to  refraction  is  computed  as 


Aa  =  it>0  -  a 


and  the  iteration  procedure  stops. 


Figure  4.5  shows  uncompensated  elevation  errors  due 
to  refraction  in  y  radians  as  a  function  of  the  actual  object 
elevation  angle  with  different  ranges  (from  500  km  to  ®)  from 
a  sensor  at  an  altitude  of  15  km.  The  atmospheric  model  used 
here  is  the  1962  U.S.  standard  atmosphere  provided  by  LOWTRAN 
5.  In  the  case  where  the  range  information  of  the  target  is 
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Fig.  4.5  The  uncompensated  elevation  errors  due  to 
atmospheric  refraction. 


not  available,  the  elevation  angle  is  compensated  for  a  preset 
range,  say  1000  km  in  the  example  shown  in  Figure  4.6.  Bias 
errors  for  targets  at  ranges  of  500,  1000  and  «  km  are  shown 
in  Figure  4.6  for  the  sensor  presented  in  Figure  4.5.  For 
other  atmospheric  models,  bias  errors  due  to  mismatched  range 
follow  similar  trends  and  magnitudes  as  those  with  the  1962 
U.S.  standard  profile;  they  are  shown  in  Figure  4.7.  For  a 
mismatched  atmospheric  profile,  Figure  4.8  shows  bias  errors 
for  objects  at  a  common  range  of  1000  km  for  all  atmospheric 
models  which  are  compensated  by  the  1962  U.S.  standard  profile 
with  the  true  range.  The  maximum  refraction  bias  that 
includes  range  and  model  mismatches  is  presented  in  Figure  4.9 
with  thereference  model  and  range  indicated.  This  figure 
presents  a  spread  of  elevation  biases  over  a  wide  range  of 
atmospheric  mismatches.  For  actual  applications,  the 
atmospheric  profile  is  known  to  a  certain  accuracy. 

Mismatches  presented  in  Figure  4.9  will  only  occur  in  a  few 
extreme  cases.  A  statistical  model  of  the  atmospheric  profile 
should  be  used  to  evaluate  the  elevation  errors  in  a 
statistical  sense. 

4.3  A  Statistical  Model  For  Atmospheric  Uncertainties 

It  would  be  very  complicated  to  express  atmospheric 
temperature  and  pressure  variations  with  a  statistical  model. 
To  evaluate  the  elevation  bias  due  to  refraction,  we  only  need 
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// 


R  =  3C 


COMPENSATED  AT  1000  km  RANGE 
WITH  MATCHED  ATMOSPHERIC 
PROFILE 


Fig.  4.7  The  bias  errors  due  to  mismatched  range  setting  in 
the  compensation  for  different  atmospheric  models. 


REFRACTION  BIAS  (jurad) 


Fig.  4.8  The  bias  errors  due  to  mismatched  atmospheric  model 
in  the  compensation  for  a  matched  range  setting. 
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MAXIMUM  REFRACTION  BIAS  (urad) 


Tfxm 


to  know  the  pressure  to  temperature  ratio.  It  is  much 
simpler  to  model  the  fluctuations  of  this  ratio.  As  we  have 
seen  in  Figure  4.3,  the  variations  of  this  ratio  are  very 
small  for  a  wide  range  of  atmospheric  models.  As  mentioned 
previously,  the  1962  U.S.  standard  atmosphere  model  has  a 
pressure  to  temperature  ratio  very  close  to  the  mean  of  the 
six  atmospheric  models  considered.  In  this  report,  we 
construct  a  statistical  model  of  the  pressure  to  temperature 
ratio  for  an  atmosphere  which  has  mean  profile  identical  to 
the  1962  U.S.  Standard  profile  as  follows: 

(P/T).  =  (P/T)  “  (z.  -z. )(A  (P/T)  +  o.o.);  i=1,2,..,34 

1  It]  1+1  1  Z  1  11 

(P/T)34  =  0,  (4.11) 

where 

iz(p/T)i  =  (PitI/Tu1  -  VV/i’i.rV 

P^  =  Atmospheric  pressure  in  millibars  of  the  1962  U.  S. 
Standard  atmosphere  at  altitude  z^ 

TV  =  Atmospheric  temperature  in  degrees  Kelvin  of  the 
1962  U.S.  Standard  atmosphere  at  altitude  z^ 
o.  =  Standard  deviation  of  A  (P/T).  which  is  a  half  of 

1  Z  1 

the  sample  standard  deviation  of  the  six  atmospher¬ 
ic  models  provided  by  LOWTRAN  5  at  altitude  z^ 

*  A  random  normal  distributed  number  with  mean  zero 
and  unit  standard  deviation. 
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Figure  4.10  and  4.11  show  the  mean  and  standard  deviation  of 
50  Monte  Carlo  runs  of  elevation  errors  after  compensation. 
The  same  conditions  are  assumed  here  as  those  in  Figures  4.6 


4.9.  It  is  interesting  to  note  that  the  mean  error  is  very 
close  to  the  bias  presented  in  Fig.  4.6  for  the  deterministic 
case  The  standard  deviation  is  less  than  200  urad  at  1"  ele 
vation  angle  for  the  worst  case. 

4.4  Geometric  Considerations 

From  the  results  presented  in  this  section,  it  is 
seen  that  refraction  errors  become  much  more  significant  at 
low  elevation  angles.  The  elevation  angle  is  a  function  of 
target  and  sensor  altitudes  and  the  distance  between  them  as 
shown  in  Figure  4.12. 
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REFRACTION  ERROR  STANDARD  DEVIATION  (yurad) 


Fig.  4.11  The  standard  deviation  of  elevation  errors  after 
compensation  for  a  statistical  atmosphere  model. 


ALTITUDE  (km) 


Fig.  4.12  The  elevation  angle  of  target  as  a  function  of 
target  altitude  and  range  to  a  sensor  at  15  km  altitude. 
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V. 


SUMMARY  AND  CONCLUSIONS 


A  formula  which  can  be  used  to  compute  the  refrac¬ 
tion  error  of  an  object  at  a  finite  range  in  a  given  atmos¬ 
phere  is  derived.  An  algorithm  which  obtains  an  approximate 
solution  of  the  formula  is  presented.  Sensitivities  of  the 
bias  error  to  the  range  and/or  atmospheric  model  used  in  the 
compensation  scheme  are  also  presented.  A  statistical  model 
which  simulates  the  effects  of  the  refraction  errors  due  to 
fluctuations  in  the  atmospheric  temperature  and  pressure  is 
presented.  Simulation  results  are  presented.  From  various 
sensitivity  studies  presented  here,  it  is  important  to 
conclude  that  the  most  important  portion  of  the  bias  error  is 
due  to  the  range  mismatch.  The  bias  error  can  be  minimized  by 
carefully  selecting  a  reference  range  within  the  operational 
range  of  the  sensor. 

A  major  conclusion  from  the  results  presented  above 
is  that  refraction  errors  are  most  serious  at  low  elevation 
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